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Bildgebungsverf ahren 

Die Erfindung betrifft ein Verfahren zur bildgebenden Unter- 
suchung eines Untersuchungsobjekts, das fur medizinische Zwe- 
cke insbesondere ein Patient ist. Das Verfahren ist insbeson- 
dere zur Anwendung im Rahmen eines Tomographie-Verf ahrens o- 
der zur Anwendung bei einer bildgebenden tomographief ahigen 
Untersuchungseinrichtung, beispielsweise bei einem Rontgen- 
computertomographiegerat, geeignet . 

Kontrastmittel sind beispielsweise bekannt aus DE 44 33 564 
Al, WO 00/16811 Oder DE 100 02 939 CI. 

Ergebnis radiographischer Verfahren, wie beispielsweise der 
Computer-Tomographie, der Mammographie, der Angiographies der 
Rontgen-Inspektionstechnik Oder vergleichbarer Verfahren, ist 
zunachst die Darstellung der Schwachung eines Rontgenstrahls 
entlang seines Weges von der Rontgenquelle zum Rontgendetek- 
tor (Projektionsbild) . Diese Schwachung wird von den durch- 
strahlten Medien bzw. Materialien entlang des Strahlengangs 
verursacht, sodass die Schwachung auch als Linienintegral tt- 
ber die Schwachungskoeffizienten aller Pixel entlang des 
Strahlenweges verstanden werden kann. Insbesondere bei Tomo- 
graphie-Verf ahren, beispielsweise bei der Rontgencomputerto- 
mographie, ist es mittels sogenannter Rekonstruktionsverf ah- 
ren moglich, von den projizierten SchwSchungsdaten auf die 
Schwachungskoeffizienten (p) der einzelnen Pixel zurttck zu- 
rechnen und damit zu einer erheblich sensitiveren Untersu- 
chung als bei reiner Auswertung der Projektionsbilder zu ge- 
langen, 

Zur Darstellung der Schwachungsverteilung wird statt des 
Schwachungskoeffizienten in der Regel ein auf den Schwa- 
chungskoeffizienten von Wasser normierter Wert, die sogenann- 
te CT-Zahl, verwendet. Diese berechnet sich aus einem aktuell 
durch Messung ermittelten Schwachungskoeffizienten ]i und dem 
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Referenz- Schwa chungskoef f izienten u Hl o nach folgender Glei- 
chung : 



mit der CT-Zahl C in der Einheit Hounsfield [HU] . Fur Wasser 
ergibt sich ein Wert Ch iO =0 HU und fur Luft ein Wert C t - - 
1000 HU. 

Da beide Darstellungen ineinander trans formierbar bzw. aqui- 
valent sind, bezeichnet im folgenden der allgemein gewahlte 
Begriff Schwachungswert Oder Schwachungskoef f izient sowohl 
den Schwachungskoef f izienten u als auch den CT-Wert. Ferner 
werden im Sachzusammenhang dieser Erf indungsbeschreibung die 
Begriff e Material und Gewebe austauschbar verwendet. Es wird 
unterstellt, dass ein Material im Kontext einer medizinisch 
angezeigten Untersuchung ein anatomisches Gewebe sein kann, 
und umgekehrt unter Gewebe in der Material- und Sicherheits- 
prtlfung ein beliebiges Material eines Untersuchungsobjekts zu 
verstehen ist. 

Obgleich die Aussagekraft eines auf den lokalen Schwachungs- 
koeff izienten (u) basierenden Bildes deutlich erh5ht ist, 
kann es im Einzelfall dennoch Probleme bei der Interpretation 
eines Bildes geben. Ein lokal erhdhter Schwachungswert lasst 
sich namlich entweder auf Materialien hoherer Ordnungszahl, 
wie beispielsweise Calcium im Skelett oder Jod in einem Kon- 
trastmittel zurtickfuhren, oder auf eine erhohte Weichteil- 
dichte, wie etwa bei einem Lungenknoten. Der lokale Schwa- 
chungskoeff izient u am Ort f ist abhangig von der in das Ge- 
webe bzw. Material eingestrahlten RSntgenenergie E und der 
lokalen Gewebe- bzw. Materialdichte p entsprechend der fol- 
genden Gleichung: 



C = 1000 x 



vy 



[HU] 



(1) 
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mit dem energie- und mater ialabhangigen Mas sen schwa chungsko- 
effizienten ^j(E,Z) und der (effektiven) Ordnungszahl Z. 

Die energieabhangige Rontgenabsorption eines Materials, wie 
sie von seiner effektiven Ordnungszahl Z bestimmt wird, tiber- 
lagert daher die von der Materialdichte p beeinflusste R6nt- 
genabsorption. Materialien bzw. Gewebe unterschiedlicher che- 
mischer wie physikalischer Zusammensetzung konnen daher im 
Rontgenbild identische Schwachungswerte aufweisen. Umgekehrt 
kann dagegen aus dem Schwachungswert einer Rontgenaufnahme 
nicht auf die Materialzusammensetzung eines Untersuchungsob- 
jekts geschlossen werden. 

Zur LSsung dieses Problems sind Verfahren zur Darstellung ma- 
terialcharakteristischer Werte erf order lich. Im Zusammenhang 
mit computerunterstiitzten Tomographieverf ahren ist es zum 
Beispiel aus US 4,247,774 bekannt, voneinander verschiedene 
Rontgenspektren oder Rontgenguantenenergien zur Erzeugung ei- 
nes Bildes zu verwenden. Derartige Verfahren werden allgemein 
als Zwei-Spektren-CT bezeichnet. Sie nutzen die ordnungszahl- 
bedingte Energieabhangigkeit des Schwachungskoef f izienten u 
aus, d.h. sie basieren auf dem Effekt, dass Materialien und 
Gewebe hoherer Ordnungszahl niederenergetische RSntgenstrah- 
lung deutlich starker absorbieren als Materialien bzw. Gewebe 
niederer Ordnungszahl. Bei hdheren Rontgenstrahlenergien 
gleichen sich dagegen die Schwachungswerte an und sind vor- 
wiegend eine Funktion der Materialdichte. Bei der Zwei- 
Spektren-CT werden dann beispielsweise die Unterschiede in 
30 den bei unterschiedlichen R6ntgenr5hrenspannungen aufgenomme- 
nen Bildern berechnet. 

Im Kontext dieser Beschreibung wird der Begriff Ordnungszahl, 
soweit nicht anders angegeben, nicht im strengen, elementbe- 
35 zogenen Sinn verwendet, sondern bezeichnet stattdessen eine 
effektive Ordnungszahl eines Gewebes, respektive Materials, 
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die sich aus den chemischen Ordnungszahlen und Atomgewichten 
der am Aufbau des Gewebes bzw. Materials beteiligten Elemente 
berechnet . 

Zu noch spezifischeren Aussagen kommt man, wenn zusatzlich 
die Methode der sog. Basismaterialzerlegung bei Rontgenauf- 
nahmen angewendet wird, wie sie etwa W. Kalender et. al be- 
schreiben in "Materialselektive Bildgebung und Dichtemessung 
mit der Zwei-Spektren-Methode, I. Grundlagen und Methodik*, 
W. Kalender, W. Bautz, D. Felsenberg, C. SUA und E.Klotz, Di- 
git. Bilddiagn. 7, 1987, 66-77, Georg Thieme Verlag. Bei die- 
sem Verfahren werden die Rontgenschwachungswerte eines Unter- 
suchungsobjekts mit RSntgenstrahlen niederer und hSherer E- 
nergie gemessen und die erhaltenen Werte mit den entsprechen- 
5 den Ref erenzwerten zweier Basismaterialien wie beispielsweise 
Calcium (fur Skelettmaterial) und Wasser (fur Weichteilgewe- 
be) verglichen. Es wird angenommen, dass sich jeder Messwert 
als lineare Superposition der Messwerte der beiden Basismate- 
rialien darstellen lasst. Zum Beispiel kann fur jedes Element 
0 der bildlichen Darstellung des Untersuchungsobjekts aus dem 
Vergleich mit den Werten der Basismaterialien ein Skelettan- 
teil und ein Weichgewebeanteil berechnet werden, so dass eine 
Transformation der ursprtinglichen Aufnahmen in Darstellungen 
der beiden Basismaterialien Skelettmaterial und Weichteilge- 
webe resultiert. 



Die Basismaterialzerlegung bzw. das Zwei-Spektren-Verf ahren 
eignen sich damit zur Auftrennung bzw. Unterscheidung von 
vordefinierten anatomischen Strukturen oder Materialarten in 
30 menschlichen und tierischen Geweben mit stark unterschiedli- 
cher Ordnungszahl . 

Aus der deutschen Patentanmeldung mit der Anmeldenummer 101 
43 131 ist ein Verfahren bekannt, dessen Sensitivitat und 
35 Aussagekraft die der Basismaterialzerlegung noch ubertrifft 
und beispielsweise eine funktionale CT-Bildgebung hoher Aus- 
sagekraft ermoglicht. Es ermSglicht die Berechnung der raum- 
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lichen Verteilung der mittleren Dichte p ( f ) und der ef f ekti- 
ven Ordnungszahl Z(r) aus einer Auswertung der spektral be- 
einflussten Messdaten einer R6ntgenapparatur . Man erhalt 
hierttber sehr gute Kontraste, insbesondere bezuglich der che- 
mischen und physikalischen Zusammensetzung des Untersuchungs- 
objekts. Beispielsweise erlaubt die Darstellung der Vertei- 
lung der Ordnungszahl im Gewebe u.a. Einblicke in die bioche- 
mische Zusammensetzung eines untersuchten Objekts, Kontraste 
aufgrund des chemischen Aufbaus in bisher dichtehomogen dar- 
gestellten Organen, eine quantitative Bestimmung von K6rper- 
bestandteilen wie z.B. Jod oder dergleichen und ein Heraus- 
segmentieren von Calcif izierungen basierend auf der Ordnungs- 
zahl . 

Es ist Aufgabe der vorliegenden Erfindung ein Verfahren an- 
zugeben, welches neue MSglichkeiten zur Empf indlichkeitsver- 
besserung oder zur ErhShung der Ausagekraft bei der mater ial- 
oder ordnungszahlabhangigen rontgenstrahlenbasierten Bildge- 
bung schafft. 

Diese Aufgabe wird gemaB der Erfindung gelost durch ein Ver- 
fahren zur bildgebenden Untersuchung eines Untersuchungsob- 
jekts, insbesondere eines Patienten, wobei 

a) dem Untersuchungsobjekt ein Kontrastmittel verabreicht 
wird, 

b) danach wenigstens zwei raumliche Verteilungen von Rontgen- 
schwachungswerten ermittelt werden, welche Rdntgenschwa- 
chungswerte jeweils den lokalen Rontgenschwachungskoef f i- 
zienten (u(x,y)) oder eine von diesem linear abhangige 
Grofie (C) darstellen, wobei die zwei raumlichen Verteilun- 
gen zumindest umfassen: 

- eine erste Schwachungswertverteilung, deren Ermittlung 
auf einem ersten Rontgenstrahlspektrum basiert, 

- eine zweite Schwachungswertverteilung, deren Ermittlung 
auf einem zweiten von dem ersten Rontgenstrahlspektrum 
verschiedenen zweiten Rontgenstrahlspektrum basiert, 
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c) unter Auswertung der beiden Schwachungswertverteilungen 
eine raumliche Verteilung eines oder mehrerer vordefinier- 
ter Ordnungszahlwerte (Z; Zl, Z2,...) oder eine raumliche 
Verteilung (Z(x,y)) nicht-vordef inierter, im Untersu- 

5 chungsobjekt vorhandener Ordnungszahlwerte ermittelt wird, 

die eine Information Uber die Verteilung des verabreichten 
Kontrastmittels im Untersuchungsobjekt enthait, und 

d) die raumliche Ordnungszahlverteilung zur bildgebenden Dar- 
stellung des Kontrastmittels verwendet wird. 

10 

Die Erfindung basiert auf dem Gedanken, dass der Einsatz von 
Kontrastmitteln die funktionale Bildgebung in der RSntgencom- 
putertomographie verbessern kann. Kontrastmittel waren hier- 
bei bislang lediglich verwendet worden, urn z.B. Blut vor sei- 
15 nem Gewebehintergrund in der Absorption abzuheben. Eine mate- 
rial- oder gewebeselektive Auswertung fand nicht statt. Der 
Erfindung liegt ferner unter anderem die Erkenntnis zu Grun- 
de, dass durch Zugabe eines Kontrastmittels in einer vertrag- 
lichen Dosis oder Konzentration eine mittels zwei verschiede- 
20 nen Rontgenspektren messbare Ordnungszahldif f erenz erzielbar 
ist . 

Bei dem Verfahren nach der Erfindung kann ein Ordnungszahl- 
wert des Kontrastmittels vordefiniert werden. Das Verfahren 
kann insbesondere mit der eingangs genannten Basismaterial- 
zerlegung kombiniert werden. 

Vorzugsweise wird die raumliche Ordnungszahlverteilung als 
zwei- oder dreidimensionales Feld ermittelt, wobei der jewei- 
30 lige Feldwert ein lokaler Ordnungszahlwert an dem durch das 
betreffende Feld reprasentierten Ort ist. Das Verfahren kann 
insbesondere mit dem Verfahren der eingangs genannten deut- 
schen Patentanmeldung 101 43 131 kombiniert werden. Der Of- 
fenbarungsgehalt dieser Patentanmeldung wird in die vorlie- 
35 gende Patentanmeldung ausdrUcklich einbezogen, insbesondere 
die dortigen Patentanspruche 1 bis 7. 
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Aufierdem bevorzugt wird zusatzlich zu der Ordnungszahlvertei- 
lung ein weiteres zwei- oder dreidimensionales Feld ermit- 
telt, dessen Feldwerte jeweils einen lokalen Dichtewert wie- 
dergeben. 

Die Verwendung der raumliche Ordnungszahlverteilung zur Bild- 
gebung kann beispielsweise dadurch geschehen, dass ein Bild 
angezeigt wird, welches nur Daten aus einem bestimmten - den 
Wert der Ordnungszahl des Kontrastmittels z.B. e ins chlieii en- 
den - Ordnungszahlintervall Oder jenseits eines bestimmten 
Ordnungszahlgrenzwertes anzeigt. Es ist auch moglich, die ge- 
messenen Ordnungszahlwerte in eine Grauwert- oder Farbskala 
umzusetzen, wobei der Wert der Ordnungszahl des Kontrastmit- 
tels ausgezeichnet oder allein eingefarbt sein kann, und die- 
se Skala bildgebend anzuzeigen. Derartigen Bildern kann ein 
gewohnliches, nicht-funktionales Schwachungsbild unter- oder 
uberlagert sein. 

Nach einer besonders bevorzugten Aus ftihrungs form werden das 
ermittelte Feld mit den Ordnungszahlwerten und das ermittelte 
Feld mit den Dichtewerten dazu verwendet, urn eine lokale Kon- 
zentration oder eine lokale Menge des Kontrastmittels zu be- 
rechnen. 

Unter einem Kontrastmittel werden im Zusammenhang mit der Er- 
findung jegliche Mittel verstanden, welche nach Zugabe in das 
Unter suchungsobj ekt, insbesondere nach Injektion in einen Pa- 
tienten, zu einer Kontrastverbesserung oder Kontrastverstar- 
kung in der Absorption, also im Rontgenbild, fuhren. Hierun- 
ter fallen sowohl konventionelle Kontrastmittel, wie sie bei- 
spielsweise bei Perfusionsmessungen in die Blutgefafie verab- 
reicht werden urn diese im Bild hervorzuheben. Es werden unter 
^Kontrastmittel' aber auch Mittel verstanden, die sich spezi- 
fisch oder selektiv, z.B. nach einem Schlussel-Schloss- 
Prinzip, nur an bestimmten Stellen im Unter suchungsobj ekt ab- 
lagern oder anreichern und somit die Uberprtifung einer Organ- 
funktion erlauben. Solche letztgenannte Mittel kSnnen auch 
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sogenannte Marker oder Tracer sein. Ein solcher Marker setzt 
sich beispielsweise aus einem biologischen Makromolektll, etwa 
einem AntikSrper, einem Peptid oder einem Zuckermolektil, mit 
einer hohen Affinitat zu der zu untersuchenden Zielstruktur, 
sowie aus einem - beispielsweise zudotierten - Kontraststof f 
zusammen, welcher im RSntgenbild gut sichtbar ist. Das Makro- 
molektll dient zum Beispiel als sogenannte "metabolische Mar- 
kierung", die bewirkt, dass sich das insgesamt auch als meta- 
bolischer Marker bezeichnete Kontrastmittel entweder aus- 
schliefilich in bestimmten Regionen, z.B. Tumoren, Entztindun- 
gen oder anderen bestimmten Krankheitsherden, anreichert. 
Kontrastmittel sind beispielsweise aus den eingangs genannten 
Schriften bekannt. 

Vorzugsweise wird ein Kontrastmittel mit einer Ordnungszahl 
grofier als 20 oder grofier als 40 verwendet wird. Das Kon- 
trastmittel weist insbesondere eine Ordnungszahl kleiner als 
83 oder kleiner als 70 auf. 

Besonders vorteilhafte Kontrastmittel enthalten Gadolinium, 
Jod, Ytterbium, Dysposium, Eisen und/ oder Wismut. 

Nach einer anderen vorteilhaften Ausgestaltung enthalt das 
Kontrastmittel eine organische Verbindung, insbesondere einen 
aliphatischen Kohlenwasserstof f , beispielsweise Zucker, und/ 
oder eine Aminosaure oder ein Peptid. 

Das Kontrastmittel kann zur selektiven Ablagerung an bestimm- 
ten Stellen oder in bestimmten Gewebeteilen des Untersu- 
chungsob j ekts ausgebildet sein. 

In vorteilhafter Ausgestaltung wird das Kontrastmittel in ei- 
ner Gewichtskonzentration aus dem Bereich 10" 4 bis 10" 7 , ins- 
besondere aus dem Bereich 10" 5 bis 10" 6 , zugegeben. 

Der im Zusammenhang dieser Schrift verwendete Begriff 'R6nt- 
genspektrum' besitzt eine weiter gefasste Bedeutung als nur 
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die Spektralverteilung einer von der Rontgenquelle der Appa- 
ratur emittierten Rbntgenstrahlung. Auch auf Seiten der Ront- 
gendetektoren kSnnen unterschiedliche Spektralanteile einer 
Strahlung mit unterschiedlichen Wirkungsgraden umgesetzt und 
somit verschieden gewichtet werden. Die daraus resultierende 
effektive Spektralverteilung wird in dieser Schrift ebenfalls 
als Rontgenspektrum bezeichnet. 

Die beiden Schwachungswertverteilungen mussen nicht notwendi- 
gerweise nacheinander als zwei Bilder mit unterschiedlicher 
Rohrenspannung auf genommen werden. Da jede RSntgenrohre ein 
Spektrum mit einer gewissen Breite emittiert, ist es bei ent- 
sprechender spektralselektiver Ausgestaltung einer zugehdri- 
gen Emp fangs einhe it auch moglich, die beiden Schwachungswert- 
verteilungen weitestgehend oder vollig simultan aufzunehmen. 
Hierzu konnten zum Beispiel in den Strahlengang zustellbare 
Filter und/ oder zwei gesonderte vorhandene RSntgendetektor- 
arrays verwendet werden. 

Insbesondere ist eine Empf angseinheit zur Durchftihrung des 
Verfahrens mit einem quantenenergieselektiven RSntgendetek- 
torarray ausgestattet . 

Insbesondere im Hinblick auf eine Verwendung des in der ein- 
gangs genannten deutschen Patentanmeldung 101 43 131 be- 
schriebenen Verfahren ist es von besonderem Vorteil, dass ei- 
ne erste funktionale Abhangigkeit eines ersten Schwachungs- 
wertes der ersten Schwachungswertverteilung von Dichte und 
Ordnungszahl und zumindest eine zweite funktionale Abhangig- 
keit eines dem ersten Schwachungswert zugeordneten zweiten 
Schwachungswertes der zweiten Schwachungswertverteilung von 
Dichte und Ordnungszahl bestimmt werden, 

und dass aus einem Vergleich der ersten funktionalen Abhan- 
gigkeit mit der zweiten funktionalen Abhangigkeit und ggf . 
weiterer funktionaler Abhangigkeiten die raumliche Ordnungs- 
zahlverteilung - und optional eine raumliche Dichteverteilung 
- ermittelt wird. 
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Vorzugsweise erfolgt hierbei die Bestimmung der funktionalen 
Abhangigkeit der Schwachungswerte von Dichte und Ordnungszahl 
ftir zumindest ein Rontgenstrahlspektrum mittels Ref erenzmes- 
5 sung an einer Eichprobe oder in Form einer Simulation auf der 
Basis eines physikalischen Modells. 

Nach einer anderen vorzugsweisen Ausgestaltung wird ein Um- 
formen der Schwachungswertverteilungen in eine Verteilung der 
10 Dichte und eine Verteilung der Ordnungszahl fUr jeden der zu- 
geordneten Schwachungswerte der ersten Schwachungswertvertei- 
lung und der weiteren Schwachungswertverteilungen auf der 
Grundlage der Ermittlung eines Wertepaares ftir Dichte und 
Ordnungszahl so vorgenommen, dass das Wertepaar die bestimm- 
15 ten funktionalen Abhangigkeiten der Rontgenabsorption von 

Dichte und Ordnungszahl far das erste Rontgenstrahlspektrum 
und zumindest ein weiteres Rontgenstrahlspektrum erftillt. Da- 
mit konnen Dichte und Ordnungszahl ftir ein Bildelement ein- 
fach als Schnittmenge der funktionalen Abhangigkeiten der 
20 einander zugeordneten Rontgenabsorptionswerte der aufgezeich- 
neten Verteilungen der Rontgenabsorptionswerte berechnet wer- 
den. 



Vorteilhafterweise weist das erste Rontgenspektrum eine Quan- 
tenenergie auf, die relativ zur Quantenenergie des zweiten 
Rontgenspektrums eine Rontgenabsorption durch den Photoeffekt 
begiinstigt, so dass eine hohe Auflosung in der Bestimmung der 
Ordnungszahlen erhalten wird. 



30 In einer bevorzugten Ausftlhrungsf orm der vorliegenden Erfin- 
dung wird zum Verandern eines Rontgenstrahlspektrums far das 
Aufzeichnen des Untersuchungsobjekts eine Veranderung zumin- 
dest eines Betriebsparameters der RontgenrShre vorgenommen, 
wobei die Rontgenquelle in einem ersten Betriebszustand ein 

35 erstes Rontgenstrahlspektrum und in einem zweiten Betriebszu- 
stand ein davon verschiedenes zweites Rontgenstrahlspektrum 
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emittiert, so dass ein schneller Wechsel zwischen zwei ROnt- 
genspektren ermSglicht ist. 

Ferner wird zum Verandern eines Rontgenstrahlspektrums fur 
das Aufzeichnen des Untersuchungsob j ekts vorteilhaft eine 
Veranderung der Detektorcharakteristik vorgenommen, wobei der 
R5ntgendetektor spektrale Teilbereiche der von der Rdntgen- 
quelle empfangenen R6ntgenstrahlung in voneinander unabhangi- 
ge elektrische Signale umsetzt und hierbei ein simultanes 
Aufzeichnen von Verteilungen der RSntgenabsorption bei unter- 
schiedlichen Rontgenspektren zulasst. 

Im Folgenden wird die vorliegende Erfindung anhand von Aus- 
fuhrungsbeispielen naher beschrieben, wobei auf folgende Fi- 
guren verwiesen wird. Es zeigen: 

Fig. 1 ein Ablaufdiagramm des Verfahrens nach der Erfin- 
dung gemafi einem Ausfuhrungsbeispiel, 

Fig. 2 anhand einer Isoabsorptionslinie das Zustandekommen 
identischer Schwachungswerte u bei Materialien un- 
terschiedlicher Zusammensetzung, 

Fig. 3a ein beispielhaf tes Funktionsschema einer Berech- 
nungsmethode zur Ermittlung von Isoabsorptionsli- 
nien als Teil des Verfahrens gemafi Figur 1, 

Fig. 3b ein beispielhaf tes Ablaufdiagramm der Transformati- 
on der Rontgenschwachungswerte in Werte der Materi- 
aldichte und Ordnungszahl als Teil des Verfahrens 
gemafi Figur 1, und 

Fig. 4 zwei Isoabsorptionslinien einer Gewebeart bei zwei 
unterschiedlichen R6ntgenstrahlspektren . 

In Figur 1 ist ein Ausfuhrungsbeispiel ftlr das Verfahren nach 
der Erfindung schematisch als Ablaufdiagramm dargestellt. In 
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einem ersten Schritt 1 wird einem Patienten P ein Tracer Oder 
ein Kontrastmittel KM verabreicht, beispielsweise durch In- 
jektion in die Blutgefafie oder durch Schlucken. Der Patient P 
wird dann in einem nur angedeuteten Rontgencomputertomogra- 
phiegerat 2 untersucht, land zwar sowohl unter Auswertung ei- 
nes ersten RSntgenspektrums SI als auch - gleichzeitig oder 
nacheinander - unter Auswertung eines zweiten RSntgenspekt- 
rums S2 (zweiter Schritt 3), die durch entsprechende Einstel- 
lung des RSntgencomputertomographiegerats 2 ausgewahlt wer- 
den. Mittels Durchftihrung einer Bildrekonstruktion (dritter 
Schritt 5) basierend auf den so erhaltenen Rohdaten wird zu 
jedem der Rontgenspektren SI, S2 eine Schwachungswertvertei- 
lung erzeugt, beispielsweise als Verteilung Ui(x,y) bzw. 
u 2 (x,y) des (linearen) Rontgenschwachungskoef f izienten u inn- 
nerhalb eines Transversalschichtbildes mit Koordinaten x und 
y. In einem vierten Schritt 7 findet computerunterstutzt eine 
Transformation der Verteilungen Ui(x,y) bzw. u 2 (x, y) des 
Rontgenschwachungskoeffizienten auf eine Ordnungszahlvertei- 
lung Z(x,y) statt. Die Ordnungszahlverteilung Z(x,y) wird in 
einem funften Schritt 8 zur Darstellung einer Verteilung des 
Kontrastmittels KM auf einem Monitor 9 genutzt. 

Geht man von einer beispielhaf ten Injektion eines Gd-basier- 
ten Tracers mit einem Gd-Atom auf 10 6 WassermolekUlen (ent- 
sprechend ca. 9 ppm Gewichtsanteil) aus, so ergibt sich ein 
effektives Z = 7,52. Verglichen mit dem Wasserwert Z = 7,42 
1st diese Konzentration mit dem Verfahren gemali der Figuren 
3a und 3b nachweisbar. Das Verfahren gemaft der Figuren 3a und 
3b kann als Teil des vierten Schritts 7 ausgefuhrt werden. 
Einzelheiten dieses Verfahrens sind in der deutschen Patent- 
anmeldung mit der Anmeldenummer 101 43 131 beschrieben, auf 
die ausfrticklich Bezug genommen wird. 

Falls gleichzeitig die Dichteverteilung p(x,y) ermittelt 
wird, kann bei bekanntem, z.B. wasserahnlichen, TrSgermateri- 
al die Konzentration des Tracerstof f es quantitativ ermittelt 
werden. 
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Die Ordnungszahl des Kontrastmittels KM sollte eine moglichst 
hohe Abweichung von der Ordnungszahl des Untergrundmaterials, 
typischerweise Wasser mit Z = 7,42, aufweisen. 

Alternativ zu dem Verfahren gemaB der Figuren 3a und 3b kann 
auch die herkc-mmliche eingangs beschriebene Basismaterialzer- 
legung verwendet werden, urn die raumliche Verteilung zweier 
oder mehr vorher festgelegter Ordnungszahlwerte Zl, Z2, . . . in 
einer Transversalschichtebene (x,y) zu ermitteln. Eine solche 
Verteilung kann ebenfalls zur Darstellung des Kontrastmittels 
KM im Bild verwendet werden. 

Im Hinblick auf eine Erlauterung des Verfahrens gemafi der Fi- 
guren 3a und 3b werden nachfolgend zunachst Erlauterungen ge- 
geben: Die Isoabsorptionslinie 11 der Figur 2 verbindet alle 
Wertepaare (p, Z) mit, bei einem definierten Rontgenstrahl- 
spektrum identischem Schwachungswert u bzw. C. Die Darstel- 
lung der Fig. 1 verdeutlicht, dass Information Uber Art und 
Zusammensetzung eines Gewebes bzw. Materials nicht allein auf 
die Schwachungswerte eines RSntgenbildes gestiitzt abgeleitet 
werden konnen. 

RSntgenstrahlung wird von unterschiedlichen Materialien und 
abhangig von der Energie der Rontgenstrahlung unterschiedlich 
stark geschwacht. Dies ist auf unterschiedlich wirkende 
Schwachungsmechanismen bei den verschiedenen Materialien zu- 
rtickzuf tihren . 

Die im Kontext dieser Beschreibung vereinfacht als Ordnungs- 
zahl titulierte effektive Ordnungszahl Z einer bestimmten Ge- 
webeart errechnet sich aus den Ordnungszahl en Zi der am Auf- 
bau beteiligten Elemente, deren Atomgewichte Ai und der en lo- 
kalen materialaquivalenten Dichten pi beispielsweise zu: 
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Z = < 



V pi 4 

Zj PiZj. 

i Al 

Z — pi 
v i Ai 



(3) 



Far reines Calcium erhalt man Z Ca =20, fUr Calciumhydrid ca. 
Z C aH2=16,04 und fur Wasser etwa Z H 2oS7, 428. Die chemische Oder 
auch biochemische Zusammensetzung eines Objekts kann daher 
sehr gut uber die Ordnungszahl Z erf ass t werden. 

Voraussetzung fur eine Berechnung der Ordnungszahl- und Dich- 
teverteilung in einem Untersuchungsgebiet sind zumindest 
zwei, in der Aufnahmegeometrie identische, aber mit unter- 
schiedlicher Energie der angewandten RSntgenstrahlung er- 
stellte RQntgenaufnahmen des Gebiets. Bei Verwendung von mehr 
als zwei mit unterschiedlicher RSntgenstrahlenergie aufge- 
zeichneten Rdntgenaufnahmen konnen die Z- und p-Aufl5sung 
verbessert werden, doch erhSht sich dadurch auch die Strah- 
lenbelastung. Im Falle der Untersuchung eines Patienten ist 
diese Moglichkeit daher nicht immer gegeben. 

Ausgangspunkt der Umwandlung von Schwachungswert basierenden 
Bilddaten in Verteilungsbilder der Ordnungszahl en und der Ma- 
terial- bzw. Gewebedichte ist die Kenntnis der Isoabsorpti- 
onslinien fur jedes RSntgenspektrum einer R5ntgenapparatur . 
Wie bereits erwahnt, ist hierbei unter ROntgenspektrum nicht 
der eng gefasste Begriff der Spektralverteilung einer von der 
RQntgenquelle der Apparatur emittierten R5ntgenstrahlung zu 
verstehen, sondern ein erweiterter Begriff, der die unter- 
schiedliche Gewichtung unterschiedlicher Spektralbereiche des 
Emissionsspektrums der RontgenrOhre auf Seiten der RSntgende- 
tektoren beriicksichtigt . Ein gemessener Schwachungswert er- 
gibt sich daher aus der direkten Schwachung des von der Ront- 
genrohre emittierten Strahlenspektrums und dem spektralen 
Wirkungsgrad des verwendeten Rontgendetektors . Beide Werte 
sind anlagenspezif ische Grofien und mUssen entweder direkt o- 
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der indirekt mittels der Schwachungswerte von Eichproben er- 
mittelt werden. Sie sind die Grundlage zur Berechnung der I- 
soabsorpt ionslinien . 

In Figur 3a sind drei Verfahren 300 zur Modellierung bzw. zur 
Berechnung einer Schar von Isoabsorptionslinien skizziert, 
namlich eine theoretische Modellierung, eine experimentelle 
Bestimmung und eine theoretische Modellierung mit einer Ka- 
librierung der Kurven durch experimentell bestimmte Parame- 
ter. 

Prinzipiell sind so viele Isoabsorptionslinien zu bestimmen, 
wie Schwachungswerte zum Abdecken der Spanne von Rontgen- 
schwachungen in den Rontgenaufnahmen erforderlich sind. Dabei 
ist nicht fur jeden theoretisch auftretenden Schwachungswert 
eine Isoabsorptionslinie zu berechnen; nicht errechnete Iso- 
absorptionslinien konnen bei Bedarf durch Interpolation oder 
andere geeignete Mittelungsverf ahren verfugbar gemacht wer- 
den. 

Die Grundschritte der theoretischen Modellierung sind im lin- 
ken Ast des Ablauf schemas der Fig. 3a dargestellt. Im Schritt 
S302 werden zunachst die Daten der fur eine Anlage spezifi- 
schen Rontgenemissionsspektren S (E) mit den verfilgbaren R5h- 
renspannungen als Parameter eingelesen. Die Spektralvertei- 
lungen der RSntgenstrahlung konnen hierzu im Vorfeld experi- 
mentell ftir jede einzelne Rontgenanlage ausgemessen werden, 
oder es werden die far einen speziellen Rontgenquellentyp 
charakteristischen Daten verwendet. Das Ermitteln der Detek- 
torapparatefunktion w(E) erfolgt in Schritt S303. Auch hierzu 
kann im Vorfeld eine genaue Vermessung der Detektoranordnung 
vorgenommen werden oder aber es werden den Detektortyp cha- 
rakterisierenden Daten wie z.B. dessen spektrale technische 
Spezifikation verwendet. Die Berechnung der Isoabsorptionsli- 
nien in Form von Kurvenscharen Ci(p, Z) bzw. Hi(p,Z) wird auf 
der Basis eines physikalischen Modells in Schritt S304 vorge- 
nommen, das fur jede relevante Kombination von S(E) und w(E) 
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die Rontgenschwachungen Ci bzw. Ui fttr Materialien mit unter- 
schiedlichen Ordnungszahlen und bei unterschiedlichen Materi- 
aldichten nachbildet. 

Alternativ zur theoretischen Modellierung der Schritte S302 
bis S304 konnen die Kurvenscharen der Isoabsorptionslinien 
auch experimentell ermittelt werden. Hierzu werden in Schritt 
S305 die Rontgenschwachungen von Eichmaterialien mit unter- 
schiedlicher Dichte und mittlerer Ordnungszahl in der Ront- 
genapparatur bei verschiedenen relevanten Kombinationen von 
S(E) und w(E) gemessen. Die Messwerte bilden die Sttitzpunkte 
ftir die folgende Berechnung der Kurvenscharen von Isoabsorp- 
tionslinien Ci bzw, Ui in Schritt S306. 

Als weitere Alternative konnen die auf theoretischer Basis 
modellierten Kurvenscharen C± bzw. ]i± mit experimentell er- 
mittelten Rdntgenschwachungswerten kalibriert werden. In 
Schritt S 307 werden die zum Eichen der theoretischen Kurven- 
scharen notwendigen Schwachungswerte wie oben ftir Schritt 
S305 beschrieben mit geeigneten Eichmaterialien bzw. Phanto- 
men in der Rontgenanlage gemessen. Im Unterschied zur rein 
theoretischen Modellierung der Schritte S302 bis S304 ist bei 
diesem Verfahren die exakte Kenntnis der Rontgenemissions- 
spektren S (E) und w(E) nicht Voraussetzung sondern Parameter 
der theoretischen Modellierung der Kurvenscharen von Isoab- 
sorptionslinien Ci bzw. ]ii in Schritt S308. Das Kalibrieren 
der Kurven in Schritt S309 mit den in Schritt S307 experimen- 
tell ermittelten Eichwerten definiert schlieJJlich Werte fttr 
diese Parameter, die spezifisch fttr die Rontgenemissions- 
spektren und Detektorapparatefunktionen der R5ntgenapparatur 
sind. 

Mit der Ermittlung der Isoabsorptionslinien fttr die erforder- 
lichen Rontgenschwachungswerte und Kombinationen von S(E) und 
w(E) sind die Voraussetzungen fttr eine Transformation von 
Bilddaten, die Schwachungswerte der Rontgenstrahlung beim 
Durchgang durch ein Gewebe reprasentieren in Bilddaten, die 
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eine Verteilung der Ordnungszahl bzw. der Mater ialdichte im 
entsprechenden Gewebe reprasentieren geschaffen. 

Je nach Auf gabenstellung kdnnen die drei Verfahren zur Isoab- 
sorptionslinienbestimmung auch gemischt verwendet werden. 
Beispielsweise konnen Werte, die experimentell nur ungenau 
oder nur mit grofiem Aufwand oder gar nicht zu ermitteln sind, 
mithilfe einer theoretischen Modellierung erganzt oder in ih- 
rer Genauigkeit prazisiert werden. Die mit unterschiedlichen 
Methoden erschlossenen Daten werden dann in Schritt S310 zu 
einem einheitlichen Datensatz zusammengef asst und in Schritt 
S311 fttr die Bildtrans format ionen bereitgehalten. 

In Figur 3b ist ein fttr das erf indungsgemafie Verfahren geeig- 
netes Transf ormationsverf ahren 320 dargestellt. Es sttitzt 
sich auf die nach einem der zuvor beschriebenen Verfahren 300 
ermittelten und als Datensatz in Schritt S321 bereitgehalte- 
nen Kurvenscharen von Isoabsorptionslinien . 

Eine Transformation erfolgt bildelementweise . Im folgenden 
wird von einer Transformation einer Rontgenschwachungswert- 
verteilung basierend auf zwei bei unterschiedlichen RSntgen- 
strahlenergiespektren aber identischer Auf nahmegeometrie auf- 
genommen RSntgenbildern ausgegangen. Dies ist die minimale 
Voraussetzung fur eine Durchftihrung einer erf indungsgemafien 
Transformation. Jedoch konnen auch mehr als zwei RSntgenauf- 
nahmen bei mehr als zwei unterschiedlichen Energieverteilun- 
gen der Rontgenstrahlung Verwendung finden. 

Die Auswahl eines zu trans formierenden Bildelements wird im 
Schritt S322 getroffen und im folgenden Schritt S323 werden 
die Schwachungswerte Ci bzw, Ui fttr dieses Bildelement aus 
dem ersten und C 2 bzw. ]i2 aus dem zwei ten R5ntgenbild gele- 
sen. Im anschliefienden Schritt S324 erfolgt die Abfrage des 
fttr die erste Rontgenaufnahme verwendeten Rontgenstrahlspekt- 
rums Si(E) und der Detektorapparatefunktionen Wi(E) sowie der 
entsprechenden Werte S 2 (E) und w 2 (E) fttr das zweite ROntgen- 
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bild. Diese Werte bilden die Parameter ftir eine nachfolgende 
Auswahl der den jeweiligen Schwachungswerten zuzuordnenden 
Isoabsorptionslinien. Die Spektralverteilungen Si(E) bzw. 
Wi(E) konnen hierbei auch indirekt, z.B. iiber eine Abfrage 
5 der verwendeten Rohrenspannungen Ui bzw. U 2 bzw. der Be- 
triebsparameter der Rontgendetektoren ermittelt werden. 

Im Schritt S325 werden aus dem in Schritt S321 bereitgehalte- 
nen Datensatz von Isoabsorptionslinien eine erste Kurve, wel- 
che die Bedingungen Ci bzw. Ui bei den Parametern Si(E) und 
Wi(E) erfullt und eine zweite Kurve, welche die Bedingungen 
C 2 bzw. u 2 bei den Parametern S 2 (E) und w 2 (E) erfullt ausge- 
wahlt. Ein Beispiel einer dergestalt erhaltenen ersten Isoab- 
sorptionslinie 11 und einer zweiten 41 Isoabsorptionslinie 
ist in Figur 4 dargestellt. 

Der Schnittpunkt 42 als Schnittmenge beider Kurven 11 und 41 
wird im Schritt S32 6 berechnet. Der Kurvenschnitt 42 lasst 
sich z.B. durch eine lokale lineare Transformation oder mit- 
tels iterativer Schnittpunktf indung ermitteln. Da die beiden 
Kurven 11 und 41 zwei unterschiedliche Schwachungswerte fur 
das selbe Bildelement und daher fur einen identischen Teilbe- 
reich eines untersuchten Gewebes reprasentieren, mtissen beide 
Schwachungswerte von der selben Material- bzw. Gewebeart ver- 
ursacht sein. Die Koordinaten (p,Z) des Kurven schnittpunkts 
42 geben daher die Materialdichte und die Ordnungszahl des 
dem Bildelement zuzuordnenden Gewebeteilbereichs wieder . 

Schliefilich wird in Schritt S327 der so ermittelte Ordnungs- 
30 zahlwert Z in die Ordnungszahlverteilung als entsprechender 
Bildelementwert geschrieben, in Schritt S328 analog der er- 
mittelte Materialdichtewert p in die Dichteverteilung. Die 
Schritte S322 bis S328 werden fur alle verbleibenden Bild- 
punkte wiederholt, bis eine abschliefiende Bildausgabe in 
35 Schritt S329 erfolgen kann. Dabei kann der Schritt S324 uber- 
sprungen werden, da die Spektralverteilungen Si(E) bzw. Wi(E) 
fur alle Bildelemente eines Bildes identisch sind. 
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Patentansprtiche 

1. Verfahren zur bildgebenden Untersuchung eines Untersu- 
chungsobjekts, insbesondere eines Patienten (P) , wobei 

a) dem Untersuchungsobj ekt ein Kontrastmittel (KM) verab- 
reicht wird, 

b) danach wenigstens zwei raumliche Verteilungen (ui(x,y), 
u 2 (x,y)) von Rontgenschwachungswerten ermittelt werden, 
welche Rontgenschwachungswerte jeweils den lokalen R6nt- 
genschwachungskoeffizienten (u(x,y)) oder eine von diesem 
linear abhangige GrOfie (C) darstellen, wobei die zwei 
raumlichen Verteilungen (ui(x,y), u 2 (x,y)) zumindest um- 

f assen: 

- eine erste Schwachungswertverteilung (Ui (x, y) ) , deren 
Ermittlung auf einem ersten Rontgenstrahlspektrum ba- 
siert, 

- eine zweite Schwachungswertverteilung (u 2 (x,y)), deren 
Ermittlung auf einem zweiten von dem ersten R6ntgen- 
strahlspektrum verschiedenen zweiten Rontgenstrahlspekt- 
rum basiert/ 

c) unter Auswertung der beiden Schwachungswertverteilungen 
(Ui(x,y), u 2 (x,y)) eine raumliche Verteilung eines oder 
mehrerer vordef inierter Ordnungszahlwerte (Z; Zl, Z2,...) 
Oder eine raumliche Verteilung (Z(x,y)) nicht-vorde- 

f inierter, im Untersuchungsobj ekt vorhandener Ordnungs- 
zahlwerte ermittelt wird, die eine Information uber die 
Verteilung des verabreichten Kontrastmittels (KM) im Un- 
tersuchungsobj ekt enthalt, und 

d) die raumliche Ordnungszahlverteilung (Z(x,y)) zur bildge- 
benden Darstellung des Kontrastmittels (KM) verwendet 
wird. 

2. Verfahren nach Anspruch 1, 

wobei ein Ordnungszahlwert des Kontrastmittels (KM) vordefi- 
niert wird. 

3. Verfahren nach Anspruch 1, 
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wobei die raumliche Ordnungszahlverteilung als zwei- oder 
dreidimensionales Feld ermittelt wird, wobei der jeweilige 
Feldwert ein lokaler Ordnungszahlwert (Z(x,y)) an dem durch 
das betreffende. Feld reprasentierten Ort (x,y) ist. 

4. Verfahren nach Anspruch 3, 

wobei zusatzlich zu der Ordnungszahlverteilung ein wei teres 
zwei- Oder dreidimensionales Feld ermittelt wird, dessen 
Feldwerte jeweils einen lokalen Dichtewert (p(x,y)) wiederge- 
ben. 

5. Verfahren nach Anspruch 4, 

wobei das ermittelte Feld mit den Ordnungszahlwerten (Z(x,y)) 
und das ermittelte Feld mit den Dichtewerten (p(x,y)) dazu 
verwendet werden, urn eine lokale Konzentration oder eine lo- 
kale Menge des Kontrastmittels zu berechnen. 

6. Verfahren nach einem der AnsprUche 1 bis 5, 

wobei ein Kontrastmittel (KM) mit einer Ordnungszahl grofter 
als 20 verwendet wird. 

7. Verfahren nach Anspruch 6, 

wobei ein Kontrastmittel (KM) mit einer Ordnungszahl groBer 
als 40 verwendet wird. 

8. Verfahren nach einem der AnsprUche 1 bis 7, 

wobei ein Kontrastmittel (KM) mit einer Ordnungszahl kleiner 
als 83, insbesondere kleiner als 70, verwendet wird. 

9. Verfahren nach einem der AnsprUche 1 bis 8, 

wobei das Kontrastmittel (KM) Gadolinium, Iod, Ytterbium, 
Dysposium, Eisen und/ oder Wismut enthalt. 

10. Verfahren nach einem der AnsprUche 1 bis 9, 

wobei das Kontrastmittel (KM) eine organische Verbindung, 
insbesondere einen aliphatischen Kohlenwasserstof f , bei- 
spielsweise Zucker, enthalt. 
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11. Verfahren nach einem der Ansprtiche 1 bis 10, 

wobei das Kontrastmittel (KM) eine Aminosaure oder ein Peptid 
enthalt . 

12. Verfahren nach einem der Ansprtiche 1 bis 11, 

wobei das Kontrastmittel (KM) zur selektiven Ablagerung an 
bestimmten Stellen oder in bestimmten Gewebeteilen des Unter- 
suchungsobjekts ausgebildet ist. 

13. Verfahren nach einem der Ansprtiche 1 bis 12, 

wobei das Kontrastmittel (KM) in einer Gewichtskonzentration 
aus dem Bereich 10 -4 bis 10"\ insbesondere aus dem Bereich 
10" 5 bis 10~ $ , zugegeben wird. 

14. Verfahren nach einem der Ansprtiche 1 bis 13, 

wobei eine erste funktionale Abhangigkeit (11) eines ersten 
Schwachungswertes der ersten Schwachungswertverteilung von 
Dichte und Ordnungszahl und zumindest eine zweite funktionale 
Abhangigkeit (41) eines dem ersten Schwa chungswert zugeordne- 
ten zweiten Schwachungswertes der zweiten Schwachungswertver- 
teilung von Dichte und Ordnungszahl bestimmt werden, 
und wobei aus einem Vergleich der ersten funktionalen Abhan- 
gigkeit (11) mit der zweiten funktionalen Abhangigkeit (41) 
und ggf . weiterer funktionaler Abhangigkeiten die raumliche 
Ordnungszahlverteilung - und optional eine raumliche Dichte- 
verteilung - ermittelt wird. 
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Zusammenfassung 
Bildgebungsverfahren 

Es ist ein Verfahren zur bildgebenden Untersuchung eines Un- 
tersuchungsobjekts, insbesondere eines Patienten (P) , be- 
schrieben. Dabei wird dem Untersuchungsobjekt zunachst ein 
Kontrastmittel (KM) verabreicht. Danach werden wenigstens 
zwei raumliche Verteilungen von RSntgenschwachungswerten er- 
mittelt, welche Rc-ntgenschwachungswerte jeweils den lokalen 
Rontgenschwachungskoeffizienten (u(x,y)) oder eine von diesem 
linear abhangige GrSfie (C) darstellen, wobei die zwei raumli- 
chen Verteilungen zumindest umfassen: 

- eine erste Schwachungswertverteilung (Ui(x,y)), deren Er- 
mittlung auf einem ersten Rontgenstrahlspektrum basiert, 

- eine zweite Schwachungswertverteilung (u 2 (X/Y)), deren Er- 
mittlung auf einem zweiten von dem ersten RSntgenstrahl- 
spektrum verschiedenen zweiten RSntgenstrahlspektrum ba- 
siert . 

Unter Auswertung der beiden Schwachungswertverteilungen wird 
dann eine raumliche Verteilung eines oder mehrerer vordefi- 
nierter Ordnungszahlwerte (Z; Zl, Z2,...) oder eine raumliche 
Verteilung (Z(x,y)) nicht-vordef inierter, im Untersuchungsob- 
jekt vorhandener Ordnungszahlwerte ermittelt, die eine Infor- 
mation uber die Verteilung des verabreichten Kontrastmittels 
(KM) im Untersuchungsobjekt enthalt. Die raumliche Ordnungs- 
zahlverteilung wird zur Darstellung des Kontrastmittels (KM) 
im Bild verwendet. 
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